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We reconstruct the shape of the primordial power spectrum of curvature perturbations in extended 
cosmological models, including addition of massive neutrinos, extra relativistic species or varying 
primordial helium abundance, from the latest cosmic microwave background data from the Wilkinson 
Microwave Anisotropy Probe, the Atacama Cosmology Telescope and the South Pole Telescope. We 
find that a scale-invariant primordial spectrum is disfavored by the data at 95% confidence level 
even in the presence of massive neutrinos, however it can lie within the 95% confidence region 
if the effective number of relativistic species or the primordial helium abundance is allowed to 
vary freely. The constraints on the extension parameters from WMAP7+ACT+_ffo+BAO, are the 
■ total mass of neutrinos < 0.48 eV (95% CL), the effective number of relativistic species 

JVoff = 4.50 ± 0.81 and the primordial helium abundance Y p — 0.303 ± 0.075. The constraints from 
WMAP7+SPT+# +BAO, are £ m v < 0.45 eV (95% CL), N cS = 3.86±0.63 and Y p = 0.277±0.050. 

^ ■ PACS numbers: 98.80.Cq 
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Measurements of temperature anisotropy in the cosmic microwave background (CMB) provide us with a wealth of 
cosmological information. Since 2001 the Wilkinson Microwave Anisotropy Probe (WMAP) satellite has measured 
| temperature anisotropy of the CMB up to I ~ 1200 over the full sky [J The ground-based telescopes such as 
the Atacama Cosmology Telescope (ACT) Q and the South Pole Telescope (SPT) have extended temperature 
Q_i measurements to smaller angular scales I ~ 3000. Measurements of the small-scale anisotropy of the CMB, in 
combination with WMAP's measurements of the degree-scale anisotropy, provide a powerful probe of early- universe 
j_i ' physics. Increasing the angular dynamic range of CMB measurements can improve the constraints on the shape of 
^ | the primordial power spectrum of curvature perturbations Moreover, hight resolution measurements of the 

small-scale CMB can limit the total mass of neutrinos, the effective number of relativistic species and the primordial 
helium abundance [!, 3 . 

^vq . Slow-roll inflationary models generically produce a nearly scale-invariant primordial power spectrum of curvature 
1 perturbations. Hence, measuring deviations from scale invariance of the spectrum is a critical test of cosmological 
OO , inflation. For a power-law parameterization of the spectrum motivated by inflationary models with featureless inflaton 
C*~) ■ potentials, the exact scale-invariant spectrum is excluded at more than 99% confidence level by the 7-year WMAP 
data 0. Other various parameterizations of the spectrum have been considered: for example, a broken spectrum Q 
caused by an interruption of the inflaton potential B, a cutoff at large scales [Tol [ill ] motivated by suppression 
of the lower multipoles in the CMB anisotropics (l2l. Il3|. and more complicated shapes of the spectrum caused by 
features in the inflaton potential [tH ]. However, a strong theory prior on the form of the primordial power spectrum 
could lead to misinterpretation and biases in parameter determination. Recently some model-independent approaches 
have been proposed to reconstruct the shape of the primordial power spectrum, based on linear interpolation 151. 
cubic spline interpolation in log-log space |g, 0] , a minimally-parametric reconstruction [loll, wavelet expansions Oj > 
principle component analysis |18| , and a direct reconstruction via dcconvolution methods In these works, the 

primordial power spectrum has been constrained within the context of the minimal ACDM model. 

On the other hand, recent CMB measurements by the ACT and SPT experiments have revealed somewhat less 
fluctuation power at small angular scales than expected in the standard ACDM model. Except for with the help of 
the primordial power spectrum with a large negative running spectral index, such a small-scale suppression in the 
CMB angular power spectrum can also be explained by extra relativistic degree of freedom or additional primordial 
helium abundance Moreover, the addition of massive neutrinos can affect the CMB anisotropies by changing 
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the matter-to-radiation ratio around the photon decoupling epoch |l| . 

In this work, we reconstruct the shape of the primordial power spectrum in the extended ACDM models, including 
addition of massive neutrinos, extra relativistic species or varying primordial helium abundance. In order to reduce 
the degeneracy between the primordial power spectrum and the extension parameters, we use the seven-year WMAP 
data in combination with small-scale CMB data from the ACT during its 2008 season and the SPT during 2008 and 
2009, and adopt two main astrophysical priors on the Hubble constant (Hq) measured from the magnitude-redshift 
relation of low-z Type la supernovae and on the distance ratios of the comoving sound horizon to the angular diameter 
distances from the Baryon Acoustic Oscillation (BAO) in the distribution of galaxies. 

This paper is organized as follows. In Section [XT] we describe the reconstruction method of the primordial power 
spectrum and the extended models in which the power spectrum is reconstructed. In Section UHl we describe the data 
used in our analysis and then present our results. Section HVl is devoted to our conclusions. 



II. BEYOND THE STANDARD MODEL 



In standard ACDM model, the primordial power spectrum is commonly parameterized as a power law, there 
are three light neutrino species in the early Universe, and big bang nucleosynthesis (BBN) took place with specific 
predictions for primordial element abundances. In this work, we consider simple extensions of the standard ACDM 
model including (i) general form of the primordial power spectrum, (ii) addition of massive neutrinos, (iii) extra 
relativistic species, (iv) varying primordial helium abundance. Each type of model is in detail described below. 

Primordial power spectrum: The primordial power spectrum, 'P(fc), is related to the angular power spectrum of the 
CMB by a radiative transfer function Af (k) via 

Cf Y = J d\nkV(k)Af(k)AY(k), (1) 

where X, Y denote the various temperature and polarization modes. Given a specific ansatz for the primordial power 
spectrum motivated by theoretical models, one can fit it to the data. Alternatively, given our complete ignorance of 
the underlying physics of the very early Universe, one can reconstruct the shape of the primordial power spectrum 
from existing data. In this paper, we use the method proposed in Q to perform a reconstruction of the primordial 
power spectrum in the extended ACDM models. We divide the spectrum into three bins equally spaced in logarithmic 
wavenumber between 0.0002 Mpc -1 and 0.2 Mpc -1 . We use a cubic spline interpolation to determine logarithmic 
value of the spectrum between knots. Outside of the wavenumber range we fix the slope of the spectrum at the 
boundaries since the CMB data place only weak constraints on them. 

Massive neutrinos: If the neutrinos are non-relativistic at the present epoch and their mass eigenstates are degen- 
erate, then the energy density parameter, f2„, can be written in terms of the total mass of neutrinos as 

n v h 2 ~ ^^r, (2) 
94eV v ; 

which as a dark matter component in the Universe cannot be negligible even for a neutrino mass as small as m v = 0.07 
cV, as established by the solar and atmospheric neutrino oscillations experiments [22|. Here, h is the dimensionless 
Hubble parameter such that Hq = lOO/i kms _1 Mpc _1 . However, at the photon decoupling epoch neutrinos were 
still relativistic as long as ^ m v < 1.8 cV, and thus the matter density at decoupling was actually smaller than a 
naive extrapolation from the present value. This leads to accelerating decay of gravitational potential around the 
decouplingepoch. Therefore, by changing the matter-to-radiation ratio the massive neutrinos modify the CMB power 
spectrum jl[ . Detecting this effect on the CMB anisotropies allows us to establish the absolute neutrino mass scale [23[ . 

Extra relativistic species: In the standard model of particle physics there are three light neutrino species. The total 
energy density of relativistic neutrinos is related to the photon energy density, p 7 , via 
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NeSP 7 , (3) 



where iV c ff is the effective number of neutrino species that is used to parameterize the relativistic degree of freedom. 
Three light neutrino species correspond to N c ff = 3.046, which slightly exceeds 3 because electron-positron annihilation 
at neutrino decoupling provides residual neutrino heating [24[ . 

The addition of extra relativistic species increases the expansion rate prior to and during the photon decoupling 
epoch and thereby increases the angular scales of photon diffusion, as explained in [25j | . Thus, increasing N e g reduces 
fluctuation power at small angular scales. The effect from extra relativistic species not only produces a nearly constant 
decrease in the angular power spectrum at I > 200 but also alters its shape at I > 200 [25j . 
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Primordial helium abundance: In the standard BBN model, the predicted value of the primordial helium abundance, 
Y p , is generally related to the baryon density and the effective number of relativistic species. For the ACDM model, 
the uncertainty in the baryon density leads to the ~ 0.02% uncertainty in the helium abundance, as N a s is fixed 
to its standard value of 3.046. Therefore, for CMB analysis the primordial helium abundance is usually assumed to 
Y p = 0.24 drl. 

Since the helium recombination happens much earlier than the hydrogen recombination, the number density of free 
electrons at the hydrogen recombination is given by n e = 715(1 — Y p ) where is the baryon number density (26|. A 
large Y p decreases the value of n e and thereby increases the mean free path of photon. This leads to lower fluctuation 
power in the damping tail of the CMB. We can extend the BBN prior by promoting Yp to a free parameter and then 
use the CMB data to probe the helium abundance. 

To summarize, we consider the extended ACDM models described by the parameters 

{n b h 2 , n c h 2 , e S! T , a u a 2 , a 3 ,x}, (4) 

where fl^h 2 and Q c h 2 are the physical baryon and cold dark matter densities relative to the critical density, Q s is the 
ratio of the sound horizon to the angular diameter distance at decoupling, t is the reionization optical depth, Ai = 
ln[10 10/ P(fci)] (i = 1,2,3) are the logarithmic values of the primordial power spectrum at knots ki, and X = ^m y , 
N e ff or Y p is the extension parameter. 

III. COSMOLOGICAL CONSTRAINTS 

In this section, we first describe the CMB data from the WMAP satellite, the ACT and SPT experiments, and 
astrophysical priors used in our analysis. We then combine these data sets to constrain the extended ACDM models 
and present our results. 

Data: We use the 7-year WMAP (WMAP 7) data including the high-/ TT power spectrum in I < 1200 and TE 
power spectrum in / < 800, and the \ow-l temperature (2 < / < 32) and polarization (2 < I < 23) data. We consider 
the Sunyaev-Zel'dovich (SZ) effect, in which CMB photons scatter off hot electrons in clusters of galaxies. Given a 
SZ template the effect is described by a SZ template amplitude Asz as in the WMAP papers [T], Q ■ 

We use the 148 GHz ACT data during its 2008 season. The focus is on using the band powers in the multiple 
range 1000 < I < 3000 to improve constraints on primary cosmological parameters. Following Q for computational 
efficiency the CMB is set to zero above I = 4000 where the contribution is subdominant, less than 5% of the total 
power. To use the ACT likelihood described in 0], aside from Asz there are two more secondary parameters, A p and 
A c . The former is the total Poisson power from radio and infrared point sources. The latter is the template amplitude 
of the clustered power from infrared point sources. We impose positivity priors on the three secondary parameters, 
use the SZ template and the clustered source template provided by the ACT likelihood package, and marginalize over 
these secondary parameters to account for SZ and point source contamination. 

We also use the 150 GHz SPT data during 2008 and 2009. We focus on the temperature power spectrum binned in 40 
band powers in the multiple range 1000 < I < 3000. Following Q Gaussian priors on the three secondary parameters 
are adopted in our analysis. The priors on the Poisson power, cluster power and SZ power are A p = 19.3 ± 3.5 /iK 2 , 
A c = 5.0 ± 2.5 /iK 2 and Asz = 5.5 ± 3.0 ^K 2 , respectively. Moreover, we require these secondary parameters to be 
positive. As discussed in [4|, constraints on cosmological parameters do not depend strongly on these priors. The 
foreground terms are used only when calculating the SPT likelihood; they are not used when calculating the WMAP 
likelihood. 

To improve constraints on cosmological parameters we follow the methodology described in [l], Q to consider the 
addition of distance measurements from astrophysical observations, on the angular diameter distances measured from 
BAO in the distribution of galaxies, and on the Hubble constant measured from the magnitude-redshift relation of 
low-z Type la supernovae. The Gaussian priors on the distance ratios are r s /Dv{z = 0.2) = 0.1905 ± 0.0061 and 
r s /Dv{z = 0.35) = 0.1097 ± 0.0036, derived from the two-degree field galaxy redshift survey and the sloan digital 
sky survey data (27j . Here r s is the comoving sound horizon size at the baryon drag epoch and Dy is the effective 
distance measure for angular diameter distance. The Gaussian prior on the Hubble constant is Hq ~ 74.2 ± 3.6 km 
s _1 Mpc^ 1 , derived from the magnitude-redshift relation of 240 low-z Type la supernovae at z < 0.1 p8j ]. The error 
includes both statistical and systematical errors. 

Results: We use a modified version of the publicly available CosmoMC package to sample the parameter space 
by means of Monte Carlo Markov Chains technique [29|. We also use the WMAP likelihood, the ACT likelihood 
and the SPT likelihood code for parameter estimation. Figure [T] shows the reconstructed primordial power spectrum 
of curvature perturbations in the standard ACDM model (top-left panel) and non-standard models with massive 
neutrinos (top-right panel), with varying relativistic species (bottom- left panel) and a varying primordial helium 
abundance (bottom-right panel), respectively, derived from WMAP7+ifo+BAO combined with the ACT small-scale 
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FIG. 1: Reconstruction of the primordial power spectrum of curvature perturbations, V(k), from WMAP7+ACT+_ffo+BAO, 
assuming — 0, N e a = 3.046 and Y p — 0.24 (top-left panel); — free, N e g = and Y v — 0.24 (top-right panel); 

]Pm„ = 0, Ngg — free and Y p = 0.24 (bottom-left panel); ^2m v = 0, N c g — 3.046 and Y p = free (bottom-right panel). 



CMB data. We can see that in the standard model the scale-invariant spectrum is disfavored by the data at 95% 
confidence level. Including massive neutrinos enhances the primordial spectrum at low-fc but suppresses it at middle- A; 
due to an enhancement in the so-called early integrated Sachs- Wolfe effect. Hence, the spectrum becomes flat at large 
scales. However, there is no significant degeneracy between massive neutrinos and the primordial power spectrum at 
small scales. Therefore, the scale-invariant spectrum is marginally disfavored by the data in the presence of massive 
neutrinos. In the extended model with free number of relativistic species (or free abundance of primordial helium), the 
scale-invariant spectrum can lie well within the 95% confidence region because the considerable degeneracy between 
relativistic species (or helium abundance) and the spectrum results in large errors and an enhancement of the spectrum 
at high-fc. Compared to varying relativistic species, varying helium abundance leads to larger uncertainties of the 
primordial spectrum at both low-fc and high-fc. Varying helium abundance affects the primordial spectrum at small 
scales by modifying the recombination process. It indirectly affects the primordial spectrum at large scales due to the 
correlation between A\ and A3 as shown in Figure O 

In Table U we list the estimated values of the primordial power spectrum as well as the extension parameters beyond 
the standard model. For a flat ACDM model with a power-law power spectrum, the WMAP7+7?o+BAO limit on the 
total mass of neutrinos is £)m„ < 0.58 eV (95% CL) 0. In our model we find that the WMAP7+ACT+ J ff +BAO 
limit is m„ < 0.48 eV (95% CL). The upper limit on the total mass of neutrinos is improved by adding the small- 
scale ACT data. With the WMAP7+ACT+i?o+BAO data combination, the effective number of neutrino species is 
estimated to be N e g = 4.50 ±0.81 (68% CL), which indicates that a Universe with no relativistic neutrinos is excluded 
at 5.6(7 and the standard model with N c n = 3.046 is consistent at 1. 8c. This estimated value is consistent with the 
effective number of relativistic species estimated in the power-law ACDM model from the WMAP7+i?o+BAO data, 
N oS = 4.34 ± 0.88 (68% CL) [|, and from the WMAP7+ACT+ J ff +BAO data, N cB = 4.56 ± 0.75 (68% CL) 0). 
Therefore, the constraints on the effective number of relativistic species do not depend strongly on additional degree 
of freedom of the power spectrum. For CMB analysis the primordial helium abundance is ususally assumed to be 
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ACDM+P(fc) ACDM+-P(fc) ACDM+P(fc) ACDM±P(fc) 

+ E?n„ + Nes + Y p 

Ai 3.054 ±0.159 3.114 ±0.165 3.039 ±0.162 3.114 ±0.172 

A 2 3.179 ±0.035 3.164 ±0.036 3.188 ±0.033 3.173 ±0.034 

A z 3.002 ±0.049 3.004 ±0.049 3.079 ± 0.060 3.049 ± 0.073 

£m„ <0.48eV 

N cS 3.046 4.50 ±0.81 3.046 

Yp (X24 (X24 0.24 0.303 ± 0.075 



TABLE I: Mean values and marginalized 68% confidence level for the primordial power spectrum parameters and extension 
parameters, derived from WMAP7±ACT±ffo+BAO. For the total mass of neutrinos, the 95% upper limit is given. 



Y p = 0.24. When it is allowed to vary, we find Y p = 0.303±0.075 (68% CL) using the WMAP7±ACT±iJ ±BAO data. 
This mean value is lower than predicted in the power-law ACDM from the WMAP7±ACT data, Y p = 0.313 ± 0.044 
(68% CL), because suppressed power spectrum of primordial curvature perturbations at small scales contributes to 
enhanced Silk damping of the CMB angular power spectrum. However, the degeneracy between the power spectrum 
and the helium abundance leads to larger uncertainties in Y p compared to the power-law spectrum. 



ACDM+P(fc) ACDlVLKP(fc) ACDM±-p(fc) ACDM±P(fc) 
+ J2 m » + N cfi + Y p 



A! 


3.048 ±0.143 


3.074 ±0.147 


3.046 ±0.147 


3.117 ±0.171 


A 2 


3.176 ±0.032 


3.168 ±0.033 


3.181 ±0.032 


3.170 ±0.034 


A, 


3.007 ± 0.040 


3.004 ± 0.040 


3.067 ± 0.059 


3.043 ± 0.063 







< 0.45 eV 








Aeff 


3.046 





3.86 ±0.63 


3.046 


Y v 


0.24 


0.24 


0.24 


0.277 ± 0.050 



TABLE II: Mean values and marginalized 68% confidence level for the primordial power spectrum parameters and extension 
parameters, derived from WMAP7±SPT±_H r o±BAO. For the total mass of neutrinos, the 95% upper limit is given. 

The constraints on the power spectrum and extension parameters from the combination data of 
WMAP7±SPT±i? ±BAO are shown in Figure H and listed in Table HH Compared to the ACT data, the addi- 
tion of the SPT data reduces uncertainties in the power spectrum at small scales. The scale-invariant spectrum is 
disfavored at 2er even in the extended model with massive neutrinos. The conclusion that the scale-invariant spectrum 
is not excluded at 2a confidence level if either the effective number of relativistic species or the helium abundance 
is allowed to vary, is confirmed by using the WMAP7±SPT±7?o±BAO data. The combination data constrain the 
total mass of neutrinos to be m v < 0-45 eV (95% CL), which is almost consistent with the result listed in Tabic Q] 
We find that the WMAP7±SPT±iJ ±BAO limit on the effective number of relativistic species is N eS = 3.86 ± 0.63 
(68% CL). In the power-law model, the constraint improves to be N g = 3.86 ± 0.42 (68% CL) jij]. Compared to the 
ACT data, the SPT data prefers to a small effective number of relativistic species. If the helium abundance is free 
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FIG. 3: Reconstruction of the primordial power spectrum of curvature perturbations, V(k), from WMAP7+SPT+7/o+BAO, 
assuming — 0, N c s — 3.046 and Y p = 0.24 (top-left panel); — free, N e g = and Y v — 0.24 (top-right panel); 

]Pm„ = 0, Ngg = free and Y p = 0.24 (bottom-left panel); ^2m v = 0, N c g — 3.046 and Y p = free (bottom-right panel). 



to vary, we find the constraint is Y p = 0.277 ± 0.050 (68% CL) from WMAP7+SPT+# +BAO. This mean value is 
lower than the result in the power-law model presented in Y p = 0.300 ± 0.030 (68% CL). Moreover, compared to 
the ACT data, the SPT data prefers to a low helium abundance. 

IV. CONCLUSIONS 

It is known that addition of massive neutrinos can influence the CMB spectrum by changing the matter-to-radiation 
ratio around the photon decoupling epoch, that extra relativistic species can suppress the CMB damping tail by 
increasing the expansion rate prior to and during the decoupling epoch, and that excessive helium abundance can also 
suppress the CMB power at small angular scales by decreasing the number density of free electrons at the hydrogen 
recombination. Therefore, the extension parameters can mimic the shape of the primordial power spectrum and 
meanwhile the general form of the primordial power spectrum can affect the constraints on the extension parameters. 

In such extended ACDM models, we have reconstructed the shape of the primordial power spectrum of curvature 
perturbations by using the cubic spline interpolation method in log-log space. To reduce the degeneracy between 
the primordial power spectrum and the extension parameters, we use the CMB data from WMAP7+ACT and 
WMAP7+SPT in combination with measurements of the Hubble constant and the BAO feature. Even if massive 
neutrinos are included in the early Universe, the scale-invariant spectrum is disfavored by the combined data at 95% 
confidence level. If the effective number of relativistic species (or the primordial helium abundance) is allowed to vary, 
then the scale-invariant spectrum can lie within the 2er bound because of the considerable degeneracy between A3 and 
N e g (or Y p ). Moreover, we found that the two data sets yield nearly the same limit on the total mass of neutrinos. 
Compared to the SPT data, the ACT data mildly prefers to a large effective number of relativistic species or a large 
primordial helium abundance. 
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The contribution from tensor modes was considered when testing for the shape of the primordial power spectrum 
in Q . Under the assumption of a scale-invariant tensor spectrum, including tensor modes can suppress the spectrum 
at large scales. In this case the reconstructed spectrum deviates from the scale-invariant one at 95% confidence level. 
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